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Abstract 



We summarize the general formalism describing surface flows in three-dimensional space in a form which is 
suitable for various astrophysical applications. We then apply the formalism to the analysis of non-radial 
perturbations of self-gravitating spherical fluid shells. 

Spherically symmetric gravitating shells (or bubbles) have been used in numerous model problems especially 
in general relativity and cosmology. A radially oscillating shell was recently suggested as a model for a variable 
cosmic object. Within Newtonian gravity we show that self-gravitating static fluid shells arc unstable with 
respect to linear non-radial perturbations. Only shells (bubbles) with a negative mass (or with a charge the 
repulsion of which is compensated by a tension) are stable. 



Subject headings: gravitation — hydrodynamics — instabilities — stars: oscillations supcrnovac: general 
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1. Introduction 



It is interesting to see how modelling problems by thin shells whose thickness is being ignored is employed in 
so many different sciences as general relativity, astrophysics, cosmology, elasticity or chemical engineering. 
In general relativity thin spherical shells of dust or perfect fluids 3 have frequently been used to analyse 
basic issues of gravitational collapse, in both its classical and quantum aspects (see e.g. Barrabes & Israel 
1991; Friedmann, Louko & Winters-Hilt 1997, and references therein); in astrophysics expanding spherical 
shells model supernovas (e.g. Vishniac 1983; Sato & Yamada 1991); the chief motivation to study shells in 
cosmology has been observations of bubble-like structures in the distribution of galaxies (e.g. Peebles 1993; 
Turok 1997) but also the physics of the early universe in which a region of false vacuum is separated by a 
domain wall (modelled usually as a spherical shell) from a region of true vacuum (e.g. Blau, Guendelman 
& Guth 1987; Bcrezin, Kuzmin & Tkachev 1987; Kolitch & Eardley 1997); in elasticity the theory of rods 
and thin shells goes back to the last century (cf. Love 1944), and in chemical engineering the mathematical 
description of the dynamics of an interface is important in such problems as calming of water waves by oil 
or in distillation and liquid extraction (Scriven 1960). A theoretical physicist of the new age would of course 
add membranes (or rather D-branes) moving in higher-dimensional spacetimes in superstring theories. 

We, as relativists, worked on various problems connected with thin shells (e.g. Bicak & Lcdvinka 1993; 
Hajfcek & Bicak 1997). Recently, one of us investigated non-radial oscillations of static self-gravitating 
spherical fluid shells in general relativity and their Newtonian limit (Schmidt 1998). However, we could 
not find a reference to this problem solved within Newtonian theory. Radial oscillations of spherical shells 
in which gravity is balanced by the surface pressure were analyzed both in the Newtonian and relativistic 
case by several authors, even for shells surrounding a compact object (Brady, Louko & Poisson 1991, and 
references therein). Most recently, in this journal such radially oscillating shells have been suggested as 
a model for variable cosmic objects (Nunez 1997). Non-radial oscillations are more difficult and even in 
Newtonian theory require some differential geometry because, for example, the correct form of the equation 
of continuity for a surface flow depends on the second fundamental form of an embedded surface in R 3 . 

In this work we investigate the non-radial oscillations of self-gravitating (or charged) spherical shells in 
Newtonian theory in detail. In Section 2 we review the formalism needed to describe surface flows. We 
here essentially follow the exposition given by Aris (1989) in chapter 10 of his book which, in turn, "is in 
the nature a somewhat extended gloss" on a paper by Scriven (1960) (both works thus emanating from a 
chemical engineering department). Our discussion is, of course, much shorter, however, it generalizes both 
mentioned works in two respects. We define a "coordinate system fixed in the surface moving in space" and 
we show how the equation of motion and the continuity equation get modified if other (general) coordinates 
are used within the surface since in concrete problems the "fixed" (Gaussian) coordinates are not practical 
at all. Secondly, although we use index notation we make occasionally contact with the formulation of 
mathematical elasticity theory by Marsden & Hughes (1983) which is based on the index-free formulation 
of modern differential geometry. In fact, Marsden and Hughes' text also touches on shell theory but it does 
not give the equations of motion and whenever it uses coordinates these are again only the " Gaussian-type" 
coordinates. 4 

In Section 3 we first discuss the general dynamics of a self-gravitating shell which is topologically spherical 
but may largely deviate from a sphere. We believe that Section 2 and the beginning of Section 3 may serve 
as a basic formalism for analyzing, for example, expanding or collapsing non-spherical sclf-gravitating shells 
in Newtonian theory in astrophysically realistic situations. In the second part of Section 3 we derive the 
conditions for a spherical shell to be in equilibrium. 

Section 4 is devoted to the derivation of the equations of motion and the continuity equation for linear 

3 By "surface perfect fluids" we mean the surface distributions of matter with isotropic stress distribution 
tangent to the surface. 

4 In Marsden & Hughes (1983), the comprehensive work by Naghdi (1972) is quoted as the standard reference 
for shells. As much as this work may be preferable for dealing with many aspects of elasticity problems, for 
our purposes, we found Aris (1989) more useful. 
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perturbations of the static solution obtained in Section 3. 

In Section 5, the stability of the static solution is analyzed. We prove that although, with an appropriate 
equation of state, the shell is stable with respect to radial oscillations, it is unstable if it is perturbed 
non-radially. It can thus hardly serve as a model for variable cosmic objects as suggested recently (Nunez 
1997). 

At the end we notice the fictitious, but amusing case of shells with negative gravitational and inertial mass. 
Since the time of an interesting work by Bondi (1957) it is well-known that in principle a negative mass can 
exist in the sense that it is not forbidden by classical physics. Some of its amusing properties were recently 
described by Price (1993). We show that spherical static shells with negative mass are, in fact, stable with 
respect to non-radial oscillations! Although there is no evidence that a negative mass exists in the real 
universe, in numerical relativity spacetimes containing negative mass solutions of Einstein's equations serve 
as testbeds. 

Finally, by considering formally the gravitational constant to be negative, we show that charged shells, in 
which the repulsive effects of the charges is compensated by a tension, are stable. We also give intuitive 
physical arguments for the results of the stability analysis in all three cases. 



2. Equations of motion for surface fluids 



The flow in a surface is more complicated than an infinitely extended 3-dimensional flow because the 2- 
dimcnsional surface (shell) can move in the 3-dimensional space which surrounds it. 

Let (t, x l ) be inertial coordinates in Newtonian spacetime. The metric gik is the time-independent metric 
on the flat Euclidian 3-space R 3 (depending on the application one may use Cartesian, polar or other 
coordinates). 

Functions x % = x l (t, a"), a = 1, 2, decribe the world lines of the particles of the fluid; a a are thus comoving 
(Lagrangian) coordinates. We assume that for fixed times t the points x l (t, a a ) form a 2-surface St in 
Euclidian space. We may think of the shell as a 3-surface in 4-dimensional spacetime which is formed by 
the flow lines, or as a sequence of 2-surfaces S t in R 3 . 

The space component of the tangent vector to the curves in 4-space, i.e. the velocity of a particle of the 
fluid in R 3 is 

' dx l 

V = I 

and its acceleration 



dt 2 



(2) 



We are free to use arbitrary coordinates z a on S t . (Objects intrinsic to S t will have Greek indices; in R 3 
we use Latin indices.) In particular we will use coordinates y a , obeying the condition that the velocities of 
the points with y"=const are orthogonal to S t for any t. In these coordinates we describe the shell by 

x* = r(t,y a ). (3) 

Then the vectors in R 3 given by 



are tangent to S t and the velocities of the points y Q =const, (df l /dt) y , in R 3 are perpendicular to t l a : 

9iA {jpk = . ( 5 ) 
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where is the metric in R 3 . These coordinates can be constructed by drawing the orthogonal curves to the 
family of 2-surfaces S t in R 3 . Choosing some coordinates y a on one surface and taking y a constant along 
the orthogonal congruence defines the coordinates y a . These coordinates are unique up to a transformation 
y a (y 13 ), independent of t. (In general relativity such a coordinate system is called "with vanishing shift" - 
see, e.g., Wald 1984.) 

In general, it may be more convenient to use other coordinates in E t , say z a , and describe the moving S t in 
R 3 by 

x* = C(t,z a ), <=(^-\, (6) 



dz a 

although {dQ /dt) z is not perpendicular to r % a — imagine, for example, a motion of the spherical surface into 
a highly oblate ellipsoidal surface which is described by 

r = R(t,0,<p) (7) 
in the standard spherical coordinates in R 3 , with z a ~ (6, ip). 

We assume y a = y a (t, a 13 ) and, inversely, a 13 = d^(t,y a )] the same for z a . The 2-dimensional metric, 
h a 0{t,y 7 ), determines the line element in S t , 

dl 2 = Kpdy a dyP , (8) 

which can be considered as induced by (pull-back of ) the metric of S t , 

dp dp 

Kf3 = dy^W 3 9l ° ' { ) 

similarly for z a . The surface velocity of the fluid is defined by 

V=(?£\, (10) 



dt 



and the acceleration by 

f dV c 



V dt 



y + V a \pV^ (11) 



where the vertical bar denotes the covariant derivative with respect to h aj 3. Such quantities defined anal- 
ogously in general surface coordinates z a do not have a natural geometrical (physical) meaning as in the 
expressions (10) and (11). 

In order to see this, imagine a particle of the fluid with fixed a a moving in R 3 according to x 1 — x l {t, a a ). 
Its velocity U l and acceleration A 1 are given by expressions (1) and (2). Using the coordinates y a , these 
may be written as 

U l = ^+t l a V a , (12) 

and 

dU l 

A^-gf + U^V* , (13) 

where t l a , V a are given by (4) and (10), and the covariant derivative is defined by f7' Q = U l a + T l - k V" t k a . Let 
n l be the unit normal to E t . Regarding equations (12) and (5), we find 

f a Ui = V a , niU z =ni-^, (14) 
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i.e., equation (12) represents the decomposition of the velocity in i? 3 into its normal and tangential parts 
with resepect to £ t . (In geometrical language, V a is the pull-back of Ui to S t .) It is easy to see that one 
can also write 

A 1 = {n^A^n 1 +t a A a , (15) 

where A a is given by expression (11) so that it represents the particle's acceleration along £ t . Now using 
general coordinates z a for which equation (5) is not satisfied, we can still write 



V° 



dz a 



for A a analogously with equation (11), and we obtain U % in the form (see eq.(6)) 



However, we now get 



U l = — + t 1 V a . 

dt a 



T l a Ui = r a gij ( — ) Z + V a , 



(16) 

(17) 
(18) 



so that V a is not a total projection of W on St (and neither its pull-back) because (dC/dt) z has also a 
non-vanishing projection onto S t . When using such coordinates we just have to remember that the surface 
velocity of the fluid is given by the whole r.h.s. of equation (18). 

Before considering the dynamics let us point out that h a fj and k a fj are geometrical quantities depending 
only on the position of the surface x 1 = C, z (t, z a ) at a given time, but independent of the particles flow. In 
fact, the set of surfaces define intrinsically the normal velocity field Uj_(t,z a ) at each point of the surface 
at a given time. The explicit expression for this field can be given in terms of our Gaussian coordinates by 
[%r] = U±n' 1 . After introducing particles their U l can be decomposed as U % = n l U± + V^ an and V a can 

be determined intrinsically by V a — (Vt an )it a - 

Now let F(t,y a ) be any function defined on S t . Denoting by St any part of S t then one can derive the 
following analogue of Reynold's transport theorem for the material (or convective) derivative of the integral 
of F over S t : 



d_ 

dt 



=/,[(f)+4 + <-» 



dS 



(19) 



where h =det(h a p) and the dot means (d/dt) y . In particular, let a(y a ,t) be the surface density of the fluid. 
Then, if the mass of any part of S t is conserved, the transport theorem (19) implies the continuity equation 



da 



+ aV c 



where da/dt — (do-/dt) y + V a dcr / dy a is the material derivative of a. 
Let us define the external curvature tensor k a p of the surface by 

*a|/J = , 

where ri 1 is the unit normal to T, t as before, and the mean curvature is defined by 5 

H = trk aP = h af3 k aP . 



(20) 



(21) 



(22) 



Our definitions agree with those in Marsden & Hughes (1983), but not with Aris (1989): b a p — —k a p, 
F(Aris) = -\H. 



A short calculation using standard geometry (see, e.g., exercise 10.41 in Aris 1989) shows that the last term 
in equation (20) can be written as 

±=HntUj, (23) 

where tP is the space velocity of a fluid particle given by equation (12). The continuity equation (20) can 
thus be rewritten in the form 

da 



-JT + <?V a la + aHrijW = , (24) 



which is the form given by Marsden & Hughes (1983) in the Theorem 5.15 and written in the Box 5.2 in 
the Gaussian coordinate system attached to S t so that equation (5) is satisfied. In general coordinates z a 
on S t , we still obtain the continuity equation in the form (24) if, instead of V a , we substitute g a,3 T^Ui 
given in equation (18); only the expression g a/3 TpUi is the component of U l parallel to S t (geometrically the 
pull-back of Ui on E t ). 

In order to derive equations of motions for the shell one starts from balancing the rate of change of momentum 
of a portion of the shell with the total force acting on it. Let the properties of the fluid be described by the 
surface stress tensor T a/3 . For example, in the case of a Newtonian surface fluid in which the viscous stress 
depends linearly on the rate of strain, the stress tensor reads (Aris 1989) 

T a P = -pg a(3 + nS a a g aSi + eE a(3po S pa , (25) 

where p is the surface pressure, the surface deformation tensor is S pa — \g ap + \ {Vp\o + Vcr\p), F a ^ \^ = 
+ ~ 9 a/3 5A/i, and k and e are the coefficients of dilatational and shear surface viscosity. We wrote 
down expression (25) just for illustration, in the following we shall consider only ideal surface fluids, i.e., 
k = e = 0, but at the moment we leave a general T a/3 . 

Let us now assume that, besides the internal pressure, there acts a surface force, F a , per unit area of the 
fluid. We require the balance of momentum in the direction of an arbitrary smooth covariantly constant 
vector field C a in the form 

f aV a C a dS = [ F a C a dS + f T afi v fj C a dl , (26) 

"* J St J St JdSt 

where 

T a dl = T al3 vpdl (27) 

is the surface stress vector acting on a linear element dl (in dSt in E t ) with a unit normal v a . Converting 
the last integral to a surface integral by Green's theorem, and using the transport theorem (19) and the 
continuity equation (20) on the left-hand side, we obtain the intrinsic (surface) equations of motion 

aA a =T af3 ]/3 + F a , (28) 

where A a is given by equation (11). (As usual, the balance of angular momentum holds since T al3 — T^ a .) 

Finally, consider the motion in i? 3 . The external force, F*, will in general have a component normal to the 
surface S t , 

F l ^(njF^rt +? a F a , (29) 

as the acceleration given by equation (15). Notice that t^F 01 , and similarly t l a T a dl (see eq. (27)), are 
just space components of the surface external force and surface stress (in geometrical language they are the 
push-forward vectors of F a and T a dl). Starting from the balance of momentum in R 3 in the direction of 
an arbitary smooth covariantly constant vector field K % analogously to equation (26) (now with quantities 
aU\F l and t 4 Q T"%), we find the complete 3-dimcnsional form of the equations of motion: 

a A 1 = {f a T a % + F l . (30) 
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Expressing t l a ^ by using equation (21) we can write them as 

a A 1 = t a T afi \p - k al3 T af) n l + F i . (31) 
In the case of a perfect fluid the equations of motion become 

r lT]i 

a— = -f a g af3 P \p + pHtf + F* , (32) 

where H is the mean curvature (22). 

It is instructive to project these equations into the directions tangent and normal to S t . The tangential part 
is given by equations (28), the normal part becomes 

a(njA>) = P H + nj F 3 , (33) 

which demonstrates how the surface pressure influences the motion in the direction normal to the surface if 
the surface is bent in R 3 so that it has non-zero external curvature. 

The equations of motion (32) and the continuity equation (24) determine the motion of the fluid. In the 
Lagrangian approach the velocity and acceleration are given by expressions (1), (2), all other quantities are 
also functions of (t, a a ), and one seeks solutions x l (t, a a ), a(t, a a ), assuming some equation of state p = p(<r) 
and external force F l are given. In the Eulcrian description one solves for f l (t, y a ) and a(t, y a ) in terms of 
which U % , t l a , V a , g a p are determined by equations (4), (9), (10), (12). Let us recall yet that we need not 
use "Gaussian-type" coordinates y a attached to the surface so that equation (5) is satisfied. We can solve 
for functions £ J (i, z a ) (see eq. (6)), we only have to remember that the surface velocity of the fluid (as it 
appears, for example, in the continuity equation (24)) is given by the whole right-hand side of equation (18). 



3. The self gravitating shell 



The gravitational field determined by a surface distribution of matter has a potential $ which is continuous 
at the surface (see, e.g., Kellog 1967). Assume that the shell S t is topologically a sphere. The derivatives of 
$ have limits on both sides of the shell, and the (covariant component of) gravitational force at a point of 
the shell is given by the mean 

Fi = , (34) 

where + $ are the potentials on the two sides of S and commas denote partial derivatives (see, e.g., 
Purcell 1965 for a derivation of equation (34) in the analogous electric case). 

Hence, if the shell of perfect fluid moves under its own gravitational field, the continuity equation has the 
form (24), and the equations of motion (32) become 

ATji 1 

= -^9 afj P\p+pHn^ - -a^i+Zj +"*, ,) . (35) 

In order that these equations indeed determine the motion of the shell we need to know $ in terms of the 
shell's variables. 

Let us assume that, in general, the position of the shell is given in spherical coordinates by r = R(t,6,ip), 
as in equation (7). Then we want to solve the Poisson equation 

Id 2 15 9$ 1 

r or* r l sin o0 o0 r 2 sin 9 dip 2 - 
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where the matter density is non vanishing only on S t . Introducing the surface matter density a(t 7 6, <p), we 
find 

P (t, r, 9, v) = a(t, 8, v) R ^^ nd S(r - R(t, 9 M ) , (37) 

where 5 is the Dirac delta and h =det(h a p) , with h a p being the metric on T, t induced by the spatial metric 
gij as in equation (9). Equation (37) follows from the general relations for 3-dimcnsional distribution Sf 
with support on a 2-surface £, given by F(r, 9, ip) = 0, which for any nice function /(r, 0, <p) requires 

/ f{r 1 9 1 i P )8 F r 2 sm6drd6dip = [ f| E rfS , (38) 
Jr 3 Jt, 

where /| E is the restriction of / to S. We shall return to the solution of equation (36) in the next section. 
Now we consider the simplest case — that of spherical symmetry. 

We thus assume a,p independent of 9,<p; U l = (R, 0,0), and the mean curvature of a sphere of radius 
r = R(t) is 

H = | , (39) 
the normal n % = (1,0,0), the potential vanishes inside and reads 

+$ = (40) 

r 

outside the shell. The equations of motion (35) thus reduce to 

aR = -ft - n^-ftr . ( 41 ) 



and the continuity equation (24) is 



R 2 R 2 



da 2<jR „ . , 

+ — = . (42) 



dt R 

The continuity equation can be immediately integrated to yield 

aR 2 = a Rl , (43) 

where the right-hand side denotes values at a fixed time (say t = 0). It is evident that equation (43) means 
the conservation of mass; it is also instructive to see how the outward directed radial force due to the surface 
pressure, i.e. the term 2p/R in equation (41), can be derived from elementary considerations of the force 
acting on a surface element dS = R 2 sin 9d9dip from its surrounding. 

Substituting M — AiraR 2 , an equation of state p = p(a), and eliminating a, by equation (43), from the 
equation of motion (41), we can solve equation (41) for the radius of the shell R(t). We shall discuss radial 
oscillations in the next section. Now we just notice that equation (41) admits a unique static solution for a 
given M and R = R. In this static case the pressure is 

_ 1 GaM 1 GM 2 

P= 4 R = = ~ R3 - ( 44 ) 

The positive surface pressure is needed to balance the inward directed gravitational force. 

It is amusing to observe that exactly the same static situation arises if M < 0, a < 0. How is it possible 
that a negative gravitational mass which repels all masses and thus pushes the elements of the shell outwards 
is compensated by a positive pressure which, as we saw above, exerts, apparently, an outward directed 
radial force on each element? The resolution of this paradox comes from the fact that we assumed both 
gravitational and inertial mass of the shell to be negative. A non-gravitational force (as the pressure) acting 
in a given direction on a negative inertial mass gives it an acceleration in the opposite direction (just by 
F = ma)! So the positive pressure, in fact, accelerates the elements of the shell in the inward direction and 
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is just compensated by the repulsive gravitational action of the negative gravitational mass. We shall see 
how these effects can influence the stability of the static shell in the following sections. 

Another alternative is to change the sign of G. Gravity becomes then electricity in which the charges of the 
same sign repel each other. Equation (44) says that we need a negative p, a tension as in soap bubbles, to 
balance such a shell. 



4. The linearized equations of motion 



In general the solution of the coupled system of nonlinear equations (24), (35) and (36) describing the 
motion of a highly deformed shell is complicated. Our goal here is to investigate linearized perturbations 
(oscillations) of a static spherical shell satisfying the conditions (44). Since the background solution is 
spherically symmetric we can, without loss of generality, assume that the perturbations are axisymmetric, 
i.e., independent of the azimuthal coordinate <p. The position S t of the shell at time t is thus given by 

r = R(t,6), (45) 

with t fixed. If we let 9 and t change, equation (45) describes a 3-surface in (Newtonian) spacetime. Now 
the fluid can move in the shell; denote its surface velocity by 

V e (t,6) = W(t,0) . (46) 

The vector 

^ = ( R, e , 1, 0) (47) 

is tangent to E t . The fluid velocity in space, given by U % = (dx l /dt)g + T l g V 6 (cf. eqs. (6 ), (17)), reads 

U l {t,6) = (R + R.eW, W, 0) , (48) 

where R = dR/dt. The surface metric ^g a /3 induced on S t is ^ggg — R 2 + R 2 , ^g vv = i? 2 sin 2 6, so that 
the covariant component V$ = (R 2 g + R 2 )W; hence, 

Tg Ui = Rfi{R + RfiW) + R 2 W , (49) 

and equation (18) is indeed satisfied. 

The acceleration appearing on the left-hand side of equation (35) is (see eq. (13)) 
which in spherical coordinates implies 

rjTjr 

— =R + 2R g W + RgW + W(RgW)e-RW 2 , (51) 
at 

^ = W + W^ + WW.g + 2W 2 ^ . (52) 
dt R R 



10 



(Since V , U e are coordinate components, their dimension and thus the dimension of W is s 1 , whereas U r 
has the usual dimension cm s _1 in CGS units.) The normal to the shell is given by 

"' = 7TTCT (1 --^' 0) - <53) 

"' - TttCT " • • 0) • <54) 

and the mean curvature turns out to be 



H = 



R^/l + {Re/Rf 



1. . 1 f R g\ f Rj 



2-^Recote + R, ee ) + i + {R9/R)2 ^) ^ 



(55) 



We can write down the exact form of the continuity equation (24) in the general case by substituting for 
rij, U 3 and H the expressions above and for V a the expression (49) since the surface coordinates 9, tp are 
generalized coordinates as z a in Section 2 rather than y a (the lines (0, <p) =const are not perpendicular to 
the surface r = R{t 1 9)). However, we shall not analyze the general case further, we shall now linearize both 
the continuity equation and the equations of motion around the static solution satisfying equation (44). 

To derive the linearized equations we consider 1 -parameter families R(t,0,e), W(t,9,e) of shell solutions. 
The coordinates in the shell and in the embedding space are uniquely fixed. Hence we obtain a description of 
the linearized equations in a particular (coordinate) gauge. As always we assume that the family is smooth 
in e and that we can interchange e-derivatives and spacetime derivatives. 

Let R(t, 9, 0) = R, W(t, 9, 0) = be a static shell. We denote background quantities with an overbar and 
the perturbation of a quantity Q by SQ. It is easy to see that in the linearized case the coordinates become 
Gaussian in linear order, i.e. as coordinates y a used in Section 2; indeed neglecting higher-order terms, 
equation (5) is satisfied. The linearized acceleration (51), (52) is 

SA l = ( 6R, W, 0) (56) 

because the background acceleration vanishes. We have now to determine the radial and tangential compo- 
nents of the linearization of all terms in equations (35) . For the inner forces due to the pressure gradient we 
obtain only a tangential component 

-Sfogofpj) = -tU^5p,p = -jpSi,8p,e (57) 

because the background pressure is independent of 9. The perturbation of the normal force is 

SipHrt) = SpHn* + pSHn 1 + pHSn 1 . (58) 
Regarding (53), (54) and (55) we obtain 

5m = ( 0, -5R,e, 0) , (59) 
Sn* = {0,-R- 2 5Rj, 0) , (60) 

SH = -j^SRfie + cot e{-^5R fi ) - -^SR . (61) 

The linearized equations of motion (35) can thus be written in the form 

ad A 1 = -8 l e R- 2 5 Pfi + Hn l 5p + pn l 6H+pHSn l 

-Isv^i+Qj - , (62) 
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where i, j — r, 8, g rr — 1, g ee = R 2 (and we do not need the (^-components because of axial symmetry). 

Before calculating the perturbations of the gravitational potential let us make the standard assumption that 
all quantities can be decomposed into spherical harmonics. Thus, we write 

oo oo 

1=0 1=0 

and 

oo oo 

SR = J2ti (t)Yi , SW = J2m (W,g , (64) 
i=o i=o 

where Yi = Yi (8), Y^g = 8Yi /d8, the form of SR and SW corresponds to the fact that 5R describes a shift 
whereas SW is a ^-component of a velocity (cf. eq. (46)). Because Y^gg + cot 9 Y^g = —1(1 + l)Yj, we obtain 
from equation (61) 



(65) 



The perturbations of the potential can be calculated by integrating the Poisson equation (36) with p obtained 
by perturbing the 5-function source (37). We shall proceed somewhat differently but we checked that both 
procedures lead to the same result. Decompose the potential inside and outside the shell into spherical 
harmonics (we are now omitting the argument t since it is irrelevant here): 

oo oo 

-$ = J2 ai(e)r l Yt , +$ = ]T MO*-'" 1 *! , (66) 

1=0 1=0 

oo oo 

-$, r - aiW^Yi , +<& r = M«0H - l)r- l - 2 Y . (67) 

1=0 1=0 

At the shell, r = R(8, e), the potential is continuous, 

-*[R(6,e),e]= + *[R(0,e),e], (68) 

and its gradient satisfies 

{n*[+<f(r, 8, e),i $(r, 8, e),,]} r=i?(M = ^Ga{8, e) . (69) 
Linearization of these relations with SR — ' ^(^' 0) = anc ^ a; = 0, fo; = for Z > 1, implies 

S ai R l = S^R- 1 - 1 - b Y R-% , (70) 

and 

-(I + l)5biR- 1 - 2 + 2b Q Y Q R- 3 £i - ISaiR 1 - 1 = inGSai , (71) 
where Y = l/V&r. We can solve for Sai, 5bi in terms of £z and 5ar. 

Sai = 2ZTT ~ ^YoR' 1 ' 2 ^ - ^GR~ l+1 S m ) , (72) 

Sk = ^— ((Z + 2)6 y i?'- 1 6 - AirGR^S^) . (73) 



For the static background shell we have b Y = — ^-kGR 2 g. Using this we obtain 

4ttC 
27+ 



Sai = WTT [ (/ " mi ~ RSai ^ R ~ l ' (74) 
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6b l = ^ I [-(l + 2)a£ l -R5a l ]R l+1 . (75) 
The gravitational intensity at the shell is 

Fi = -l[ + *Me), i +-*(rAe),i\ r=me) ■ (76) 

Inserting ± $ we obtain after linearisation, using that the background is spherically symmetric, for the 
covariant ^-component (with fixed angular behavior Yj) of the force intensity at the shell 

F w = -\ [S ai R l + SkR- 1 - 1 ] Y lfi , (77) 

and the radial component is 

F lr = - X - [SarfR 1 - 1 - 8bi(l + 1)R- 1 - 2 + 2b Y R- 3 Ci] Y x . (78) 

Inserting 8ai,5bi and b in which for I = wc omit the background term, —GM/2R 2 , since it drops out as 
a consequence of equation (44), we obtain 

Fw = [-3*6 - 2R5m] Y Xfi , (79) 

and 

+4irGaR- 1 &Y l . (80) 

It remains to consider the linearized form of the continuity equation (24). Since a =const, and V a — W = 0, 
the first term, da/dt = da/dt + V a da/dy a , after linearisation just becomes 8a = dSa/dt, the second term 
becomes a(8W.g + SWcot9), and the third an r 8RH. Substituting the angular decompositions (63), (64), 
and using Yi_gg + cot 9 Y^g = —1(1 + l)Yz, we obtain the Z-part of the continuity equation in the form 

8m - 1(1 + l)afn + 2aR- 1 ii = . (81) 

Integrating we get 

Sai - 1(1 + l)arn + 2aR~^i = , (82) 
where we put the integration constants equal to zero since = r\i = implies 8<ji = 0. 

Multiplying now the intensity components (79), (80) by a, assuming the background conditions (44) satisfied, 
and substituting the perturbed quantities as given above into the equations of motion (62), we obtain 

aii = 2R- 1 8 P i + jGaMR~ 3 [l(l + 1) - 2]& - \pMRT 2 8a l 

-^iflim-^R-^i + Scn], (83) 

am = -R- 2 8 Pl - ^GaMR-% + la^^R- 2 (3a^ + 2R8a l ) , (84) 
where the first equation is meaningful for all I > 0, the second for I > 1. 

Finally, let us assume that the perturbed pressure and matter densities are connected by a linear relation 
Sp = a8a, so that 

Spi = a8ai . (85) 
13 



This, by p = p(<r), 8p = (dp/ da)5a, corresponds to a general equation of state for barotropic fluids. Substi- 
tuting for Sp into (83), (84), writing a — M/AttR 2 , and expressing Sai from the continuity equation (68), we 
arrive at a system of two equations for just £j and r\i as follows: 



GM 



I 2 -2,1 -2 2-1(1 + 1) 4aR 



21 + 1 



GM 



m 



GM 



+■ 



GM 

w 



1(1 + 1) 2 



21 + 1 



I + 1 2aR 



21 + 1 GM 



+ 21(1 + 1) 



xR 



GM 



6 + 



GM 

w 



1(1 + 1) 

21 + 1 



1(1 + 1) 



aR 



GM 



m 



(86) 



(87) 



where for I = only the first equation is meaningful. Since we already used the continuity equation the last 
two equations are the only equations to be solved for , m to determine the general axisymmetric linearized 
perturbations. 



5. The stability analysis 



Before investigating stability let us cast the last coupled equations (86), (87) into a still simpler form. 
Denoting 

li = 6 , m=Rm, = «^ . (88) 



and the dimensionlcss time coordinate 



we get 



r=[^) t, (89) 



dT 2 
dT 2 



Al + Bfn=0, (90) 
C£i + Df,i=0, (91) 



where the coefficients are 



. I 2 -2.1 -2 2-1(1 + 1) ao 
B=^-2l(l + l)0, 

D = - l -^+l(l + l)f3. (92) 

Applying d 2 /dr 2 to equation (90) and regarding equation (91), we obtain the following 4-th order equation 
for|/, 

^. + (A + D)^ + (AD-BC)i l =0, (93) 
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and the same equation for fji. Assuming = ^\e %UlT , the last equation implies 

ujf - (A + D)ujf + AD - BC = . 
Hence, the frequencies of the oscillations are given by 



(94) 



,(1,2) 



± r-(A+ D )± 



j(A + Df - {AD - BC) 



(95) 



where (1,2) refers to the ± sign inside the bracket; the first sign (outside the bracket) represents a trivial 
alternative. The system is stable if — given a fixed (3 — the w;'s are real for all I. 

Let us first look at radial oscillations. With / = 0we have A = — | + 4/3, B = 0, C = 1 — 2/3, D = 0, so that 
the expression (95) gives 

u l=0 = ±2^3-^J . (96) 

The second solution of the relation (94), ui — 0, has no meaning since for I = only the first equation (92) 
with t) — is valid. Therefore, we conclude that the shell is stable with respect to radial oscillations if 

»\. 

Using equation of state 5p = a Sa and relations (44), (88), we can write this as the condition 



a 



Sp 3GM 
> 8R 



a 5v 3 
<=> — — > — 
pSa > 2- 



(97) 



This displays the analogy to the standard stability condition, Ti = (p/p)(Sp/Sp) > 4/3, for radial adiabatic 
stellar oscillations. For stronger gravity (3GM/8R large), stiffer equation of state, Sp — a 5a, is needed to 
guarantee stability. 

Turning next to dipole (/ = 1) perturbations we get A = — | +4/3 = — B, C = | — 2/3 = — D, and the formula 
(95) yields 







,(2) 



2(3/3- 1) 1/2 . 



(98) 



Regarding the equations (90), (91), we easily find out that the second solution implies a trivial amplitude 
|i = fji = 0, whereas the first, ojj^ = 0, implies time independent amplitudes; by incorporating the angular 
parts we easily see that, as usually, they just describe a (small) shift of the origin of the coordinates along 
the axis 6 = 0, it. 

For quadrupole (/ = 2) perturbations we get 



-(A + D) 



+ 5/3 



2 V ' 2 

so that the equation for the frequencies, (95), implies 



AD-BC 



(99) 



J l=2 



(-^ + 5/3) 2 + ^ 



(100) 



,(2) 



=2-± -2 +5/? ~ 



(101) 



Since for stable radial oscillations we must have /3 > 



(2) 



is not real and therefore the self-gravitating 



fluid shell is unstable with respect to quadrupole pertubations. 
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Although this result is of course sufficient to prove instability let us look at perturbations with large I. We 
find equations (92) to imply for I — > oo 

\(A + D)^ l -l\(3-\)+ l -l((3+\) + 0(l) , 



AD-BC 



-l 2 (3 + 0{l) . 



(102) 



Again, /3 > | implies \{A + D) > but AD — BC is negative so that uj®^ is imaginary — there is an 
instability. 

Let us now turn to the amusing case of shells with negative mass M (and with also a < as the inertial 
mass desity). Putting first I = in the original equation (86) for £/, and defining (3 again by (88), we obtain 



GM 

w 



4/3 Uo = . 



(103) 



Introducing dimensionless time 



we find 



-GM 



LU l=0 = ±2 ^-/3 + 
Since /3 = aR/GM, we get stability for any (3 such that 

3G(-M) 



-oo < 13 < 



i.e. 



8R 



< a < oo . 



(104) 



(105) 



(106) 



Introducing r from equation (104) into the expressions (86), (87), we obtain again the equations (90), (91), 
only with negative signs of A, B, C, D given by (92). Thus, the frequencies of oscillations of the shells with 
M < are 



(1,2) 



±\-\{A + D)± 



-{A + Df - {AD-BC) 



(107) 



The difference between this expression and expression (95), i.e. the sign change at the first term ^(A + D), 
is crucial. For dipole perturbations nothing new arises — one solution for uj leads to vanishing amplitudes, 
the other represents just a shift of the origin. However, turning to the quadrupole oscillations we now get 
frequencies 



<!-««•+!* 



which are all real provided that 



P < , i.e. a > . 

Similarly, regarding the expressions (102) for large I, and taking f3 < we find all frequencies real: 



,(!) 



±1/(1-4/3)5 , 



(2) 



~ ±- 



/3 



"1-4/3 



(108) 
(109) 

(110) 



In fact, one can prove that for any integer I > 1 the expression in superbrackets under the square root in 
equation (107) is positive. Indeed, the expression under the square root in the square brackets is a parabola 
in (3 for any fixed I and one can show that its minimum is always positive. 
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Next, one finds that 



l(A + D) = - 1 -(l 2 + l + 4)p+ 1 -(f + l + 6) 



which is positive for (3 < 0. Hence the expression with the plus sign in the superbrackets is positive. By 
forming its product with the corresponding expression with the minus sign one obtains 



AD - BC = 



4(2/ + 1) 1 



-/3(2Z-3)(Z + 2) + Z-2] 



(1 2) 

which is always positive for integer I > 1 (though not e.g. for I = 1.2), and so uj\ ' ' is always real. 

Thus, we conclude that self-gravitating spherical shells with negative both inertial and gravitational mass are 
stable with respect to small oscillations. 

Finally, it is of interest to consider an analogous case in electrostatics. We thus neglect all inductive and 
radiative properties of the electromagnetic field, taking only into account the fact that the same charges repel 
each other by Coulomb's law. This should be a reasonable approximation for charged bubbles. Hence, we 
put M — Q > 0, the total charge of the sphere, and consider a negative gravitational constant — G = 7 > 0. 
In fact, we can take any negative value for 7, multiply by 7 _1 both original equations (83) and (84), and 
identify then by 

£=*, = - (ill) 

7 



7 



the inertial surface mass density and its perturbation, whereas leaving a — Q/AttR 2 and 5a as the background 
charge density and its perturbation. (We thus consider a fluid of charged particles with a fixed value of the 
specific charge.) 

First, from the equilibrium condition (43) we find that in order to have a static equilibrium we need a 
negative pressure, i.e. a tension 

Q 2 



p = -7 



14i? 4 



which balances the repulsive electric force. Next we assume again a linear relationship 

5pi = ctSai = ajSHi . 



(112) 



(113) 



The equation of continuity, (82), is valid for both inertial mass and charge. However, in the terms on the 
r.h.s. of the equations of motion (83), (84), giving the electric force, we of course have to substitute 5ai. It 
is easy to see that the equations (86) and (87) become 



7^ = - 



Q_ 



Q_ 

R 2 



1.. Q 



1 



l 2 -31 -2 
21 + 1 

1(1 + I) 2 
21 + 1 

2aR 



1(1 + 1) AaR 



Therefore, writing 



21 + 1 7 Q 
III, fji = Rqi as in equation 



P Q 



and defining now 



m ■ 



(114) 
115) 



aR aR 



(116) 



and 



iQ 
i? 3 



t , 



(117) 
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we arrive at 



dr 2 



Mi 



Bfji 



, 



(118) 




C6 



Dm 



o , 



(119) 



dr 2 



where A, B, C, D are given again by equation (92). These are exactly the same equations as those we analyzed 
for the shells with negative mass. Consequently, we can conclude that the charged shells will be stable with 
respect to radial and non-radial oscillations if the parameter (3 in equation (116) is negative. Since G < 0, 
this requires 



i.e., regarding equation (113), a decrease of the tension when the mass and the charge density increase. 
Consider a static charged shell in which the repulsive effects of the charges is compensated by a tension. If 
the radius of the shell is increased while its total charge is unchanged, the repulsion decreases and a decreased 
tension can still bring back the shell into the original radius. 

The results of the stability analysis of all three problems we discussed can, in fact, be made intuitively 
plausible by considering the following Figure 1. 

In Figure la a perturbed self-gravitating shell with positive mass is illustrated. Gravity is pointing always 
inwards. But the perturbed element is pushed inwards also by the surface pressure exerted by adjacent 
elements, so an instability arises. In Figure lb a shell with a negative mass is considered. Gravity now 
acts outwards, the pressure acts also outwards but it causes an acceleration pointing inwards because the 
inertial mass of the element is negative (see the discussion at the end of Section 3), and thus the shell with 
negative mass can be in stable equilibrium. If the element is pushed inwards as in Figure la, both gravity 
and pressure give it an acceleration pointing outwards. With charged shells the situation is similar to that 
of a shell with a negative mass. It is now a tension which pulls an element inwards in Figure lb and thus 
acts against the electric repulsion. 

In the theory of adiabatic nonradial stellar oscilations the Schwarzschild discriminant, A = (dhip/dr) — 
T^ 1 (dlnp/dr), plays an important role. If the criterion of convective stability, A < 0, is violated in the 
whole star, unstable modes exist (e.g. Ledoux & Walraven 1958). This corresponds to our considerations 
illustrated in Figure 1. In the situation described in Figure la the element suffers a "convective instability", 
whereas it is stable in the situation depicted in Figure lb. 
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Figure 1 

a) A perturbed self-gravitating shell with positive mass. Both gravity and pressure on the perturbed element 
point inwards. The shell is unstable. 

b) A perturbed self-gravitating shell with negative gravitational and inertial mass. Both gravity and pressure 
on the perturbed element point outwards but the acceleration caused by the pressure points inwards due to 
the negative inertial mass. A similar situation arises with charged shells with tension. Such shells are stable. 
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